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Abstract 

A classification of discrete integrable systems on quad-graphs, i.e. 
on surface cell decompositions with quadrilateral faces, is given. The 
notion of integrability laid in the basis of the classification is the three- 
dimensional consistency. This property yields, among other features, 
the existence of the discrete zero curvature with a spectral parameter. 
For all integrable systems of the obtained exhaustive list, the so called 
three-leg forms are found. This establishes Lagrangian and symplectic 
structures for these systems, and the connection to discrete systems of 
the Toda type on arbitrary graphs. Generalizations of these ideas to 
the three-dimensional integrable systems and to the quantum context 
are also discussed. 



1 Introduction 

The idea of consistency (or compatibility) is in the core of the integrable 
systems theory. One is faced with it already at the very definition of the 
complete integrability of a Hamiltonian flow in the Liouville-Arnold sense, 
which means exactly that the flow may be included into a complete family 
of commuting (compatible) Hamiltonian flows Similarly, it is a charac- 
teristic feature of soliton (integrable) partial differential equations that they 
appear not separately but are always organized in hierarchies of commuting 
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(compatible) systems. It is impossible to list all applications or reincarna- 
tions of this idea. We mention only a couple of them relevant for our present 
account. A condition of existence of a number of commuting systems may 
be taken as the basis of the symmetry approach which is used to single out 



integrable systems in some general classes and to classify them ||2J]. With 



the help of the Miwa transformation one can encode a hierarchy of integrable 



equations, like the KP one, into a single discrete (difference) equation ^5 
Another way of relating continuous and discrete systems, connected with the 
idea of compatibility, is based on the notion of Bdcklund transformations 
and the Bianchi permutability theorem for them |Q. This notion, born in 
the classical differential geometry, found its place in the modern theory of 
discrete integrable systems A sort of the synthesis of the analytic and 
the geometric approach was achieved in and is being actively developed 
since then, see a review in These studies have revealed the fundamen- 
tal importance of discrete integrable systems; it is nowadays a commonly 
accepted idea that they may be regarded as the cornerstone of the whole 
theory of integrable systems. For instance, one believes that both the differ- 
ential geometry of "integrable" classes of surfaces and their transformation 
theory may be systematically derived from the multidimensional lattices of 
consistent Backlund transformations [Q. 

So, the consistency of discrete equations steps to the front stage of the 
integrability theater. We say that 

a (i-dimensional discrete equation possesses the consistency 
property, if it may be imposed in a consistent way on all d- 
dimensional sublattices of a (d + l)-dimensional lattice 

(a more precise definition will be formulated below). As it is seen from the 
above remarks, the idea that this notion is closely related to integrability, 
is not new. In the case d = 1 it was used as a possible definition of inte- 
grability of maps in [Bl]. A clear formulation in the case d = 2 was given 



recently in [^. A decisive step was made in [Q: it was shown there that 
the integrability in the usual sense of the soliton theory (as existence of the 
zero curvature representation) follows for two-dimensional systems from the 
three-dimensional consistency. So, the latter property may be considered 
as a definition of integrability, or its criterion which may be checked in a 
completely algorithmic manner starting with no more information than the 
equation itself. Moreover, in case when this criterion gives a positive re- 
sult, it delivers also the transition matrices participating in the discrete zero 
curvature representation. (Independently, this was found in PB| .) 

In the present paper we give a further application of the consistency 
property: we show that it provides an effective tool for finding and classi- 
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fying all integrable systems in certain classes of equations. We study here 
integrable one-field equations on quad-graphs, i.e. on cellular decomposi- 
tions of surfaces with all faces (2-cells) being quadrilateral In Sect, ^we 
formulate our main result (Theorem 1), consisting of a complete classification 
of integrable systems on quad-graphs. Of course, we provide also a detailed 
discussion of the assumptions of Theorem 1. Sect. 0, H are devoted to the 
proof of the main theorem. In Sect. | we discuss the so called three-leg forms 
PI of all integrable equations from Theorem 1. This device, artificial from 
the first glance, proves to be extremely useful in several respects. First, it 
allows to establish a link with the discrete Toda type equations introduced in 
full generality in . Second, it provides a mean to the most natural deriva- 
tion of invariant symplectic structures for evolution problems generated by 
equations on quad-graphs. This is discussed in Sect. ^. Further, Sect. ^ 
contains a brief discussion of the relation of the equations listed in Theorem 

1 with Backlund transformations for integrable partial differential equations. 
Finally, in Sect. |^ we discuss some perspectives for the further work, in par- 
ticular the consistency approach for three-dimensional systems, and discrete 
quantum systems. 

2 Formulation of the problem; results 

Basic building blocks of systems on quad-graphs are equations on quadrilat- 
erals of the type 



Q{x,u,v,y;a,p) = 0, 



(1) 



where the fields x,u,v,y & C are assigned to the four vertices of the quadri- 
lateral, and the parameters a,(3 E C are assigned to its edges, as shown on 
Fig.§ 




Figure 1: An elementary quadrilateral; fileds are assigned to vertices 
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A typical example is the so called "cross-ratio equation" 



{x-u){y-v) ^ a 
{u-y){v-x) /?' 

where on the left hand side one recognizes the cross-ratio of the four complex 
points X, u, y, v. We shall use the cross-ratio equation to illustrate various 
notions and claims in this introduction. 

Roughly speaking, the goal of the present paper is to classify equations 
([l|) building integrable systems on quad-graphs. We now list more precisely 
the assumptions under which we solve this problem. 

First of all, we assume that the equations (|l]) can be uniquely solved for 
any one of its arguments x,u,v,y G C. Therefore, the solutions have to 
be fractional-linear in each of their arguments. This naturally leads to the 
following condition. 

1 ) Linearity. The function Q{x, u, v, y; a, (3) is linear in each argument 
(affine linear): 

M, f , y; a, (3) = aixuvy + ■ ■ ■ + aig, (3) 

where coefficients depend on a, (3. 

Notice that for the cross-ratio equation (0) one can take the function on 
the left-hand side of (|l|) as m, f , y; a, [3) = P{x—u){y—v)—a{u—y){v—x). 

Second, we are interested in equations on quad-graphs of arbitrary com- 
binatorics, hence it will be natural to assume that all variables involved in 
the equations (|I]) are on equal footing. Therefore, our next assumption reads 
as follows. 

2) Symmetry. The equation (|lD is invariant under the group 1^4 of the 
square symmetries, that is function Q satisfies the symmetry properties 

Q{x, u, V, y; a, (3) = eQ{x, v, u, y; (3, a) = aQ{u, x, y, v; a, (3) (4) 

with e,a = ±1. 

Of course, due to the symmetries (^) not all coefficients in are 
independent, cf. formulae (123) below. 



We are interested in integrable equations of the type (|l|), i.e. those ad- 
mitting a discrete zero curvature representation. We refer the reader 
to 0, 0, §1, where this notion was defined for systems on arbitrary graphs. 
As pointed out above, in the third of these papers it was shown that the in- 
tegrability can be detected in an algorithmic manner starting with no more 
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Figure 2: D4 symmetry 



information than the equation itself: the criterion of integrabihty of an equa- 
tion is its three-dimensional consistency. This property means that the 
equation (|I|) may be consistently embedded in a three-dimensional lattice, 
so that similar equations hold for all six faces of any elementary cube, as on 
Fig. H (it is supposed that the values of the parameters aj assigned to the 
opposite edges of any face are equal to one another, so that, for instance, all 
edges (a;2,a;i2), (0:3, Xsi), and (X23, 3:123) carry the label ai): 




Figure 3: Three-dimensional consistency 



To describe more precisely what is meant under the three-dimensional 
consistency, consider the Cauchy problem with the initial data x, Xi, 0:2, X3. 
The equations 

Q{x ) = (5) 

allow one to determine uniquely the values X12, 0:23, X31. After that one has 
three different equations for 0:123, coming from the faces (xi, X12, xsi, X123), 
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(a;2, a;23, 3^12, 3^123)5 and (xa, X31, X23, X123). Consistency means that all three 
values thus obtained for X123 coincide. 

For instance, consider the cross-ratio equation (^). It is not difficult to 
check that it possesses the property of the three-dimensional consistency, and 

^ _ (ai - ^2)3^13^2 + (^3 - ai)x3Xi + {a2 - ^3)3^23^3 
(^3 - «2)3;i + («! - a3)3;2 + («2 " ai)3;3 

Looking ahead, we mention a very amazing and unexpected feature of the 
expression (||): the value 2:123 actually depends on a;i,X2,3;3 only, and does 
not depend on x. In other words, four black points on Fig.|^ (the vertices of 
a tetrahedron) are related by a well-defined equation. One can rewrite the 
equation (||) as 

(xi - X3){x2 - X123) _ ai - a3 



(xa - X2) (3:123 - 3;i) 02 -«3 



which also has an appearance of the cross-ratio equation for the four points 
(xi, 0:3, X2, 3:123) with the parameters ai — assigned to the edges (^1,^3), 
(0:2, 3:123) and a2 — as assigned to the edges (0:2,3:3), {xi, 3:123). 

This property, being very strange from the first glance, holds actually not 
only in this but in all known nontrivial examples. We take it as an additional 
assumption in our solution of the classification problem. 

3) Tetrahedron property. The function X123 = z{x, Xi, X2, x^; ai, 02, as) 
existing due to the three-dimensional consistency, actually does not depend 
on the variable x, that is, = 0. 

Under the tetrahedron condition we can paint the vertices of the cube into 
black and white ones, as on Fig. ^, and the vertices of each of two tetrahedrons 
satisfy an equation of the form 



Q{Xi, X2, Xs, X123; «!, "2, "3) = 



it is easy to see that under the assumption 2) (linearity) the function Q may 
be also taken linear in each argument. (Clearly, formulas (|^), (^ may be 
also written in such a form.) 

Actually, the tetrahedron condition is closely related to another property 
of Eq. (|ip, namely to the existence of a three-leg form of this equation 

V^(x, u; a) - iIj{x, v; (3) = 0(x, y; a, (3). (9) 

The three terms in this equation correspond to three "legs" : two short ones, 
{x,u) and {x,v), and a long one, {x,y). The short legs are the edges of the 
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original quad-graph, while the long one is not. (We say that the three- leg 
form (P) is centered at the vertex x; of course, due to the symmetries of the 
function Q there have to exist also similar formulas centered at other three 
vertices involved.) For instance, the cross-ratio equation (|^) is equivalent to 
the following one: 

^ ^^^^ 



u — X V — X y — X 

The three-leg form gives an explanation for the equation for the "black" 
tetrahedron from Fig.|^. Consider three faces adjacent to the vertex X123 on 
this figure, namely the quadrilaterals 0:12, X31, X123), (a;2, ^23, a;i2, a;i23), 
and (xs, 0:31, X23, X123). A summation of three-leg forms (centered at X123) of 
equations corresponding to these three faces leads to the equation 

0(xi23, Xi, a2, as) + 0(xi23, X2; as, «i) + 4>{xi23, X3; «!, ^2) = 0. 

This equation, in any event, relates the fields in the "black" vertices of the 
cube only, i.e. has the form of For example, for the cross-ratio equation 
this formula reads: 



Xl — X123 X2 — X123 X3 — X123 

and a simple calculation convinces that this is equivalent to (^). 

So, the tetrahedron property is a necessary condition for the existence 
of a three-leg form. On the other hand, a verification of the tetrahedron 
property is much more straightforward than finding the three-leg form, since 
the latter contains two d priori unknown functions ip, (j). 

It remains to mention that, as demonstrated in [^], the existence of the 
three-leg form allows one to immediately establish a relation to discrete 
systems of the Toda type Indeed, if x is a common vertex of n 
adjacent quadrilaterals faces (x, x^, x^^fc+i, Xfc+i), k = 1,2,... ,n, with the 
parameters assigned to the edges (x, x^) (cf. Fig.^), then the fields in the 
point X and in the "black" vertices of the adjacent faces satisfy the following 
equation: 

n 

^0(x,Xfc,A;+i,afe,afc+i) = 0. (11) 

k=l 

This is a discrete Toda type equation (equation on stars) on the graph whose 
vertices are the "black" vertices of the original quad-graph, and whose edges 
are the diagonals of the faces of the original quad-graph connecting the 
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"black" vertices. The parameters a are naturally assigned to the corners 
of the faces of the "black" subgraph. See Fig. ^. (Of course, a similar Toda 
type equation holds also for the "white" subgraph.) 

For instance, in the case of the cross-ratio equation, the discrete Toda 
type system (0) reads as 

n 

y^^^^ = o. (12) 

In the next Section we will show that the tetrahedron condition naturally 
separates one of two subcases of the general problem of classification of inte- 
grable equations on quad-graphs. The second subcase will not be considered 
in this paper; the corresponding subclass of equations is certainly not empty, 
but we are aware only of trivial (linearizable) examples. 

By solving the classification problem we identify equations related by 
certain natural transformations. First, acting simultaneously on all variables 
X by one and the same Mobius transformation does not violate our three 
assumptions. Second, the same holds for the simultaneous point change of 
all parameters a ^ '^ip)- Our results on the classification of integrable 
equations on quad-graphs are given by the following statement. 

Theorem 1. U'p to common Mobius transformations of the variables x and 
point transformations of the parameters a, the three- dimensionally consistent 
quad-graph equations with the properties 1),2),3) (linearity, symmetry, 
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tetrahedron property) are exhausted by the following three lists Q, H, A (u — 
xi, V = X2, y ^ Xi2, ai, (3 ^ a2). 
List Q: 

(Ql) a{x -v){u-y)- (3{x - u)(v - y) + 5'^a(3(a - (3) ^0, 

(Q2) a{x — v){u — y) — [3{x — u){v — y) + a[3{a — [3){x -\- u -\- v + y) 

-aP{a - /3)(q;2 - aP + P^) = 0, 



(Q3) (/S^ - a'^){xy + uv) + l3{a'^ - l){xu + vy) - a{l3'^ - l){xv + uy) 

-(52(o;2 - /32)(q,2 _ l)(/32 _ 1)/(4q;/3) = Q, 



(Q4) ttQXUvy + ai{xuv + uvy + vyx + yxu) + a2{xy + mv) + 02 (xm + 

+a2(a:'i; + uy) + 03(3; + rt + v + + 04 = 0, 

where the coefficients ai are expressed through (a, a) and {(3, h) with — 
r^a), = r{P), r{x) — Ax^ — g2X — g^, by the following formulae: 

ao — a + b, ai = —[3a — ab, 02 = (3'^a + a^b, 

93 92 92 

List H: 

(HI) {x -y){u-v) + p - a^O, 

(H2) (x - y)(u - v) + ((3 - a)(x + u + V + y) + (3"^ - a"^ ^ 0, 
(H3) a{xu + vy) - j3{xv + uy) + 5{a^ - p^) = 0. 
List A: 

(AI) q;(x + f;)(M + I/) - p{x + u){v + y) - S'^ap{a - /?) = 0, 

(A2) {p^ - a'^) {xuvy + 1) + p{a'^ - l){xv + uy) - a{p^ - 1) {xu + vy) = 0. 
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Remarks 1) The list A can be dropped down by allowing an extended 
group of Mobius transformations, which act on the variables x, y differently 
than on u,v (white and black sublattices on Figs.|I|,|^). In this manner Eq. 
(Al) is related in to (Ql) (by the change u — > —u, v —v), and Eq. (A2) 
is related to (Q3) with 6 = (by the change m— So, really 
independent equations are given by the lists Q and H. 

2) In both lists Q, H the last equations are the most general ones. This 
means that Eqs. (Q1)-(Q3) and (HI), (H2) may be obtained from (Q4) 
and (H3), respectively, by certain degenerations and/or limit procedures. 
So, one could be tempted to shorten down these lists to one item each. 
However, on the one hand, these limit procedures are outside of our group 
of admissible (Mobius) transformations, and on the other hand, in many 
situations the "degenerate" equations (Q1)-(Q3) and (HI), (H2) are of an 
interest for themselves. This resembles the situation with the list of six 



Painleve equations and the coalescences connecting them, cf. |jT5 . 

3) Parameter S in Eqs. (Ql), (Q3), (H3) can be scaled away, so that one 
can assume without loss of generality that 6 = or 6 = 1. 

4) It is natural to set in Eq. (Q4) (a, a) = {p{A),p'{A)) and, similarly, 
{P,b) = {p{B) , p' (B)) . So, this equation is actually parametrized by two 
points of the elliptic curve //^ = r(A). The appearance of an elliptic curve 
in our classification problem is by no means obvious from the beginning, its 
origin will become clear later, in the course of the proof. For the cases of 
r with multiple roots, when the elliptic curve degenerates into a rational 
one, Eq. (Q4) degenerates to one of the previous equations of the list Q; for 
example, if (72 = fi's = then the inversion a; — > 1/x turns (Q4) into (Q2). 

Bibliographical remarks. It is difficult to track down the origin of the 
equations listed in Theorem 1. Probably, the oldest ones are (HI) and 
(H3)^^Q, which can be found in the work of Hirota (of course not on 
general quad-graphs but only on the standard square lattice with the labels 
a constant in each of the two lattice directions; similar remarks apply also 
to other references in this paragraph). Eqs. (Ql) and (Q3)^^q go back to 



30| , see also a review in |2^. Eq. (Q4) was found in p|. A Lax represen- 



tation for (Q4) was found in |^ with the help of the method based on the 
three-dimensional consistency, identical with the method introduced in 
Eqs. (Q2) and (Q3)^^^ are particular cases of (Q4), but seem to have not 
appeared explicitly in the literature. The same holds for (H2) and (H3)^^^. 
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3 Classification: analysis 



In principle, the three-dimensional consistency turns, under the assumptions 
1), 2), into some system of functional equations for the coefficients Oj of the 
functions Q (see (0)). However, this system is difficult to analyze and we will 
take a different route. 

For the first step, we consider the problem of the three-dimensional con- 
sistency in the following general setting: find triples of functions /i, /2, /a of 
three arguments such that if 

2^23 = /l(a;,a;2,X3), Xgi = /2(a;, Xg, Xi), Xi2 = f'i{x,Xi,X2) (13) 

then 

a;i23 := /l(a;i,Xi2,X3i) = f2{x2,X2Z,Xi2) = f3{x3,X3i,X23) (14) 

identically in x,Xi,X2,X3. In other words, we ignore for a moment the con- 
ditions 1) and 2) and look to what consequences the tetrahedron condition 
leads. The proof of the following statement demonstrates that this condition 
separates just one of two possible subcases in the general problem. 

Proposition 2. For the functions /i,/2,/3 compatible in the sense 
l ^m ) and satisfying the tetrahedron condition, the following relation holds: 

f3,X2f2,xifl,X3 = — f3,xxf2,xzfl,X2 (15) 

identically in x, Xi,X2,X3. 

Proof. Denote Xus = z{x,Xi,X2,X3), and also denote for short the functions 
with shifted arguments (|T^) by capitals: Fi = fi{xi,xi2,X3i), etc. Then 
differentiating (|T^) with respect to X2,X3 and x yields the following system 
which is linear with respect to the derivatives of Fi: 

f'i,X2-Fl,Xi2 ~ ^X2 

f2,xsFi^X3i = 
f3,xFl,x\2 ~^ f2,xFl,X3i ^x 

and two analogous systems for F2, F3 obtained by the cyclic shift of indices. 
The above system is overdetermined, and, excluding the derivatives of Fi, 
we come to the first equation of the following system (the other two come 
from the similar considerations with F2, F^): 

f2,xsf3,xZx2 ~^ f3,X2f2,xZx3 /3,X2/2,X3^x 
flyXifSjX^xi ~l~/3,xi /l,a;^X3 fl,X3f3,xiZx 

fl,X2f2,x^xi ~^f2,xifl,xZx2 f2,xifl,X2^x 

Now, the tetrahedron condition = implies that the r.h.s. vanishes, and 
therefore the determinant has to vanish as well; this is equivalent to (|TB[) . □ 
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The second possibility, which we do not consider here, would be that 
^ ^ and Eq. ([TSl) does not hold. In this case the above system can be 
solved to give z^.^ = ^iZ^, % = 1,2,3, where $j are expressed through /j j,, 
fj^Xk- The compatibility of the latter three equations can be expressed as a 
system of differential equations for <l>j, and therefore for which are much 



more complicated than ([T5|). It certainly deserves a further investigation, see 
discussion in Section H. 

The necessary condition (|T5|) will provide us with a finite list of candidates 
for the three-dimensional consistency. First of all, we rewrite the relation 



(]T3|) in terms of the polynomial Q using both the linearity and the symmetry 



assumptions. 

Proposition 3. Relation ^TB) is equivalent to 



g{x, xi] ai, a2)g{x, x-i, ot^, ot^)g{x, X3; 03, ai) 

= -g{x, xx\ ai, OL'^gix, X2; ^2, ai)g{x, x^, 03, 02), (16) 

where g{x, u; a, (3) is a biquadratic polynomial in x, u defined by either of the 
formulas 

g{x,u;a,f3) = QQyv - QyQv, (17) 
g{x,v;p,a) = QQyu- QyQu, (18) 

where Q = Q{x, u, v, y; a, (3) . The polynomial g is symmetric: 

g{x,u;a,/3) = g{u,x;a,p). (19) 

Proof. The equivalence of the definitions (p!7D, (p!8| ) follows from the first 
symmetry property in (^ {a ^ (3, u ^ v), while the second one {x ^ u, 



y v) implies (19). 

To prove (0), let Q = p{x,u,v)y + q{x,u,v), so that y = f{x,u,v) = 
-q/p. Then f^/fu = {qvP - qPv)/{quP - qPu), and substituting p = Qy, 
q = Q — yQy, we obtain 

fv QQyv QyQv 



fu QQyu QyQu 

which yields: 



fk,xj _ g{ai,aj;x,Xi) 



(^,j,A;)=c.p.(l,2,3). 



fk,Xi fl'(cij; Clji Xj) 

Now ([TB|) follows from (|I§). □ 
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The biquadratic polynomials (0) and (p!8| ) are associated to the edges of 
the basic square. One can consider also the polynomial 

G{x,y;a,(3) = QQuv - QuQv, (20) 

associated to the diagonal. They have the following important property. 

Lemma 4. The discriminants of the polynomials g = g{x,u] a, (3) , g = 
g{x,v; (3,a), and G = G{x,y; a, (3) , considered as quadratic polynomials in 
u, V, and y, respectively, coincide: 

gl - 2gguu = t - ^gg^ = Gl - 2GGyy. (21) 

Proof. This follows solely from the fact that the function Q is linear in each 
argument. Indeed, calculate 

9u '^99uu i,i,QQyv QyQv^u) ^i^QQyy QyQv){,QQyv QyQv)uui 

taking into account that Q^u = 0. The result reads: 

9u '^99uu — Q Quvy ~^ QuQvy ~^ QvQuy ~^ QyQuv ~^ ^QQuvQuyQvy 
^QQuvyi^QuQvy ~t~ QvQuy ~t~ QyQuv^ ^QuQvQyQuvy 

2QuQvQuyQvy ^QnC^yC^y^-uC^yy yQ iifjQ ^^y . 

It remains to notice that this expression is symmetric with respect to all 
indices. □ 

In the formula (|16]) the variables are highly separated, and it can be 
effectively analyzed further on. In the next statement we demonstrate that 
this functional equation relating values of g with different arguments implies 
some properties for a single polynomial g. 

Proposition 5. The polynomial g{x,u; a, (3) can be represented as 

g{x, u; a, (3) = k{a, P)h{x, u; a), (22) 

where the factor k is antisymmetric, 

k{p,a) = -k{a,p), (23) 

and the coefficients of the polynomial h{x,u;a) depend on parameter a in 
such a way that its discriminant 

r(x) = hl- 2hh^^ (24) 

does not depend on a. 
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Proof. Relation (|T6|) implies that the fraction g{x, xi, ai, a2)/g{x, xi, ai, a^) 
does not depend on xi, and due to the symmetry (p^Of ) it does not depend on 
X as well. Therefore 

g{x,xi;ai,a2) K{ai,a2) 



5f(x,xi; ai, as) K{ai,as)' 
where, because of (|16D , the function n satisfies the equation 

K{ai, a2)K{a2, a3)K{as, ai) = -K{a2, ai)K{as, a2)K{ai, a^). 
This equation is equivalent to 

K{f3,a} = —^K{a,l3), 

that is, the function = is antisymmetric. We have: 

g{x,u]a,[3) _ g{x,u;a,'y) ^ g{x,u;a,(3) _ g{x,u;a,'y) 
K{a,P) K(a,7) k{a,P) k{a,'y) ' 

which implies (0). To prove the last statement of the proposition, we notice 
that, according to (^2]) , (p3D , 

k{a, P)h{x, u; a) = g{x, u; a, P), —k{a, P)h{x, v; [3) = g{x, v; P, a). 

Due to the identity (pTD, we find: 

/i^ — 2hhuu = hl — 2hhyy, h = h{x, u;a), h = h{x, v; /?), 

and therefore r does not depend on a. □ 

Thus, the three-dimensional consistency with the tetrahedron property 
implies the following remarkable property of the function Q which will be 
called property (R): 

(R) the determinant g = QQvy — QvQy is factorizable as in (^2]), (p3D, 
and the discriminant r = h^^ — "ihhuu of the corresponding quadratic 
polynomial h does not depend on parameters at all. 

It remains to classify all functions Q with this property. This will be done 
in the next section. A finite list of functions Q with the property (R) con- 
sists, therefore, of candidates for the three-dimensional consistency. The final 
check is straightforward, and shows that the property (R) is not only nec- 
essary but also almost sufficient for the three-dimensional consistency with 
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the tetrahedron property (the hst of functions with the property (R) consists 
of a dozen of items; only in two of them one finds functions violating the 
consistency). 

As a preliminary step, we consider more closely the coefficients of the 
polynomial Q. The symmetry property @) easily implies that two cases are 
possible, with a = 1 and a = — 1, respectively: 

Q = a^xuvy + ai{xuv + uvy + vyx + yxu) 

+ a2{xy + uv) + 02 (xM + vy) + 0,2 (xf + uy) 

+ a3(x + M + + ?/) + a4, (25) 

Q = ai{xuv + uvy — vyx — yxu) 

+ a2{xy — uv) + a^{x — u — v + y), (26) 

where 

ai{P,a) = eai{a, (3), 02(0, /3) = £a2(/5, a), £ = ±1. (27) 

It is easy to prove that the case (|26| ) is actually empty. Indeed, a light 
calculation shows that in this case 

g{x, u; a, (3) = g{x, u; /3, a) = —P{x; a, P)P{u] a, P), 

where P{a, /3; x) = aix"^ — a2X — 03, so that the relation ([TBI) becomes 

P{xi; ai, a2)P{x2] 02, a3)P(a;3; 03, ai) 

= -P{xi] as, ai)P{x2; ai, a2)P{x3; 02, "s)- 

Equating to zero the coefficients at {xiX2Xs)'^, X1X2X3, and the free term 
yields ai = 02 = 03 = - a contradiction. 

Turning to the case (|25|), we denote the coefficients of the polynomials 
h, r as follows: 

h(x, u; a) = hox^u^ + hixu{x + u) + 62(2;^ + u^) + h2xu + 63(0; + u) + 64, 

(28) 

r{x) = Cqx'^ + Cix^ + C2X^ + C3X + C4, (29) 

where hi = bi{a). So, we consequently descended from the polynomial Q 
(4 variables, 7 coefficients depending on 2 parameters) to the polynomial h 
(2 variables, 6 coefficients depending on 1 parameter; also take into account 
the factor k), and then to the polynomial r (1 variable, 5 coefficients, no 
parameters) . 

It remains to go the way back, i.e. to reconstruct k, h and Q from a given 
polynomial r (which does not depend on the parameter a). This will be done 
in the next section. 
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4 Classification: synthesis 



First of all, wc factor out the action of the simultaneous Mobius transforma- 
tions of the variables x. The action x ^ {ax + h)/ {cx + d) transforms the 
polynomials /i, r as follows: 



Using such transformations one can bring the polynomial r into one of the 
following canonical forms, depending on the distribution of its roots: 

• r{x) — 0; 

• r{x) — 1 (r has one quadruple zero); 

• r{x) = X (r has one simple zero and one triple zero); 

• r{x) — x^ (r has two double zeroes); 

• r{x) = x^ — I (r has two simple zeroes and one double zero); 

• r{x) = ix^ — g2X — gs, A = g(| — 27 ^ (r has four simple zeroes). 

Next, we find for these canonical polynomials r all possible polynomials h. 

Proposition 6. For a given polynomial r{x) of the fourth degree, in one 
of the canonical forms above, the symmetric biquadratic polynomials h{x, u) 
having r{x) as their discriminants, r{x) — h\ — 2hhuu, are exhausted by the 
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following list: 

r = : h = —(x — uf; (qO) 

a 

/i = (70XM + 7i(x + m) + 72)^; (hi) 

r = l: h = —{x-u) --■ (ql) 

h = '^q{x + uf + -ii{x + u) + -i2, 7i- 47072 = 1; (h2) 

r = X : h = —ix — uf — —{x + u) + '^-, (q2) 

4a 2 4 



r = : h = 7oX^u^ + 71XM + 72, 7i - 47072 = 1; (h3) 

2 r2 . " / 2 , 2^ 1 + "^ , X2l~"^ ^ ON 

r = a; — : = -{x + u ) -xu + ; [qi) 

[xu + a(x + u) + (72/4)^ 



4x^ - 5(2X - 5(3 



— (x + M + Q;)(4axM — (73) . (q4) 



Proof. We have to solve the system of the form 

hi - 46062 = Co, 

261(62 -262) -46063 =ci, 
hi - Ahj - 26163 - 46064 = C2, 

263(62 -262) -46164 =C3, 
6^ - 46264 = C4, 



where hk are the coefficients of h{x,u) and are the coefficients of r(x) (see 
(|28|), (p9|)). This is done by a straightforward analysis. For example, consider 
in detail the case r = 4x^ — g2X — g^. We have: co = 0, ci = 4, hence 60 7^ 0. 
Set 60 = 61 = — 2ap~^, then 62 = Next, use the second and the 

third equations to eliminate 63 and 64, then the last two equations give an 
expression for 62 and the constraint = r{a). □ 

Remarks. 1) Notice that for all six canonical forms of r we have a one- 
parameter family of /I's, denoted in the list of Proposition ^ by (q0)-(q4) 
(the one-parameter family for r = is not explicitly written down since 
it coincides with (q3) at 5 = 0). It will turn out that the polynomials Q 
reconstructed from these /I's belong to the list Q (so that h in (qk) corre- 
sponds to Q in (Qk); h from (qO) corresponds to (Ql)5=o)- poly- 
nomials r = 0,1, we have additional two- or three-parameter families of 
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/i's denoted in the list of Proposition ^ by (hl)-(h3). They will lead to the 
correspondent Q's in (H1)-(H3), as well as to (Al), (A2). 

2) The expression ^/r{a) in Eq. (q4) clearly shows that the elhptic curve 
/i^ = r(A) comes into play at this point. One can uniformize ^/^(a) = p'{A), 
a = p{A), where A is a point of the elliptic curve, and p is the Weierstrass 
elliptic function. Actually, the polynomial h{x, u; a) is well-known in the the- 
ory of elliptic functions and represents the addition theorem for p-function. 
Namely, the equality h{p{X), p{U); p{A)) = is equivalent to A = ±X ± U 
modulo the lattice of periods of the function p. To resume: symmetric bi- 
quadratic polynomials h{x,u) with the discriminant — 2hhuu = t{x) are 
parametrized (in the non-degenerate case) by a point of the elliptic curve 
/i^ = r{X). This is the origin of the elliptic curve in the parametrization 
of Eq. (Q4). As a consequence, the spectral parameter in the discrete zero 
curvature representation of the latter equation also lives on the elliptic curve. 
This may be considered also as the ultimate reason for the spectral param- 
eter of the Krichever-Novikov equation to live on an elliptic curve (|^, see 
also Sect.^. 

It remains to reconstruct polynomials Q for all /I's from Proposition ^ 
Proof of Theorem 1. For the polynomial (^) we have: 

g{x, u; a, (3) = (a2ao — al)x'^u'^ + (ai(a2 — 0.2) + oocts — aia2)xu{x + u) 
+ (ctiaa — 020.2) (x^ + u^) + (0,2 — a\ + aQa^ — a\)xu 
+ (03(02 — 02) + O1O4 — 0203) (a^ + u) + 0204 — 03, 

and an analogous expression for g{x, u; P, a) is obtained by the replacement 
02 ^ 5,2- Using (|2^) and denoting bi = bi{a), b[ = bi{(3), k = k{a,P), we 
come to the following system for the unknown quantities a^: 



— 2 

O2O0 — 0,1 




= kbo, 


~ 2 

O2O0 — O-]^ 




= 'kb',, 


01(02 - O2) + O0O3 - 


O1O2 


= kbi, 


01(5,2 - 02) + O0O3 - 


O1O2 


= -kb[, 


O1O3 — O2O2 




= kb2 


O1O3 — O2O2 




= -kb'2, 


O2 — O2 ~t~ O0O4 — O2 




= kb2, 


O2 — O2 ~l~ O0O4 — O2 




= -kb'2. 


03(02 — O2) + O1O4 — 


0203 


= kbs, 


03(02 — 02) + O1O4 — 


O2O3 


= -kb's. 


— 2 
O2O4 — O3 




= kbi, 


2 

O2O4 — O3 




= -kb'^. 



Of course, also the quantity k is still unknown here. Since we are looking for 
the function Q up to arbitrary factor, it is convenient to denote 



o,- '^ 02 1 02 1 k , , 

02-^2, Ai = ^, A2 = --- = - + -, K = - 30 

o 2 o 2 
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(it is easy to see that a ^ since otherwise h = 0). These functions are 
skew-symmetric: 

and the above system can be rewritten, after some elementary transforma- 
tions, as follows: 

K[{k + b'2)ibo + b'o) - (&i + b[r] = 2{bo - b',), 

K[{bo + b',){bs + b',) - (61 + 6;)(62 + b',)] = 61 - b[, 



K[{bi + b[){b, + b',) - (62 + b',){b2 + y,)] = 2(62 - b'. 



2; 5 



K[{bo + b',){b, + b',) - (62 + b',f] = ^(62 - b',), 

K[ib, + b[){b, + b',) - (62 + b',)ib, + b',)] = 63 - &3, 
K[ih + b'2)ib4 + b',) - (63 + b',f] = 2(64 - 61), 

Ao = K{bo + b'o), 2Ai = K{bi + fo'i), A2 = K{b2 + 6'2), 
4I2 = K{b2 + b'2), 2A3 = K{b^ + 6'3), = K{b^ + 

The first six lines here form a system of functional equations for bi and K 
(call it (if, 6)-system), while the last two line split away, and should be 
considered just as definitions of Ai. So, to any solution of the {K, 6)-system 
there corresponds a function Q whose coefficients are given by the last two 
lines of the system above, and the formulas 

^2 = ^2 + 1/2, I2 = I2 - 1/2, 



which follow from (30). By construction, this function has the property (R) 
and is, therefore, a candidate for the three-dimensional consistency with the 
tetrahedron property. 

Consider first the cases (qO), (ql), (q2), (q3), (q4), i.e. when we have a 
one-parameter family of polynomials h. A straightforward, although tedious, 
check proves that in these cases all six equations of the {K, 6)-system lead to 
one and the same function K, provided the functions bi = bi{a) are defined 
as in (q0)-(q4). Calculating the corresponding coefficients Ai, we come to 
the functions Q given in Theorem 1 by the formulas (Ql)5=0' (Ql)<5=i; (Q2), 
(Q3), and (Q4), respectively. A further straightforward check convinces that 
all these functions indeed pass the three-dimensional consistency test with 
the tetrahedron propery. 

It remains to consider the cases (hi), (h2), (h3). A thorough analysis of 
the {K, 6)-system shows that in these cases we have the following solutions. 

In the case (hi): 
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• either h does not depend on parameters at all, and then it has to be of 
the form h — {{^EqX + ei){eQU + £i)) , and K is arbitrary; performing 

a suitable Mobius transformation, we can achieve that h = 1\ the 
outcome of this subcase is the following function Q which is a candidate 
for the three-dimensional consistency: 

Q^{x-y){u-v) + k{a,(5), (HI) 

where k is an arbitrary skew- symmetric function, k{a.j3) = ~k{j3,a)] 

• or h = {l/a){eoXU + ei{x + u) + e2)^ with arbitrary constants Si. In this 
case we can use Mobius transformations to achieve h = {l/a){x + uy, 
and the correspondent function Q coincides with (AI)^^q from Theorem 
1. This function passes the three-dimensional consistency test with the 
tetrahedron property. 

In the case (h2): 

• either h = x + u + a, and the correspondent function Q coincides with 
(H2) from Theorem 1; 

• or h = {l/2a){x + m)^ — (a/2), and the correspondent function (Q) is 
given in (Al)^^^ of Theorem 1. 

Both functions are three-dimensionally consistent with the tetrahedron prop- 
erty 

Finally, in the case (h3): 

• cither h does not depend on parameters at all, then it has to be equal 
to h = xu, and K is arbitrary; we have in this subcase the following 
function Q which is a candidate for the three-dimensional consistency: 

1 + k 1 — k ~ 
Q = ^^{xu + vy) —{xv + uy), (H3o) 

where k is an arbitrary skew-symmetric function, k{a,(3) — —k{P,a); 

• or h — xu + a (possibly, upon application of the inversion x ^ 1/x, 

I— > l/u); the correspondent function Q is (H3)^^j^ of Theorem 1; 

a 1 + 

• or, finally, h = 7t(x'^u'^+1) rxw; the correspondent function 

1 — 1 — 

Q is given in (A2) of Theorem 1. 
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In the two last subcases the three-dimensional consistency condition is ful- 
filled with the tetrahedron property. 

To finish the proof of Theorem 1, we have to consider the functions Q 
given by the formulas (HI) and (H3o). They depend on an arbitrary skew- 
symmetric function k, and have the property (R) for any choice of the latter. 
As it turns out, these are the only situations when the property (R) does not 
automatically imply the three-dimensional consistency. 

A direct calculation shows that the function (HI) gives a map which is 
three-dimensionally consistent, with 

^ _ k{ai, a2)xiX2 k{a2, a^)x2X2, fe(Q;3, ai)x2,Xi 
k{az, a2)xx -\- k{ax, az)x2 -\- k{a2, ai)x3 

if and only if 

k{ai, Q!2) + k{a2, as) + k{a3, ai) = 0. (32) 
To solve this functional equation, differentiate it with respect to ai and 0:2: 

kaiaiioil, Oi2) = 0, 

which together with the skew-symmetry of k yields k{ai, 0:2) = f{a2)—f{oLi). 
A point transformation of the parameter /(a) ^ a allows us to take simply 
k{a,f3) = j3 — a. Thus, we arrive at the case (HI) of Theorem 1. 
Finally, consider the formula (HSq). Denote, for brevity, 

£(a,m = i±-fl^, so that £(/3, a) = l/£(a, /3). 
1 — k{a, p) 

We will write also lij for -^(aj, aj). A straightforward inspection shows that 
the function (H3o) gives a three-dimensionally consistent map with 

(^23 - 42)a;i + (4i - ^3)3:2 + (^12 - ^21)2:3 ' ^ ' 

if and only if 

£(«!, a2)^(«2, tt3)^("3, "1) = 1- (34) 

Just as above, up to a point transformation of the parameter, the solution 
of this functional equation is given by i{a, = a/P, which leads to the 
equation (H3)^^q of Theorem 1. 

The proof of Theorem 1 is now complete. □ 
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5 Three-leg forms 



Theorem 1 classifies the three-dimensionally consistent equations under the 
tetrahedron condition. In Section |l| we have seen that the latter condition is 
necessary for the existence of a three- leg form 

ipi^x, u; a) — i/j^x, v; (3) = (f){x, y; a, (3) (35) 

of Eq. (0). Recall that the formula (|35|) implies the validity of equations 
on stars ([Tl|), or equations of the discrete Toda type, for the fields x at the 
vertices of the "black" sublattice; the same holds for the "white" sublattice. 
We prove now that a three-leg form indeed exists for all equations listed in 
Theorem 1. After that, in Sect. ^ we demostrate some further applications 
of the three-leg forms, establishing the variational and symplectic structures 
for Eqs. (|) and (|Tl]). 

The theorem below provides three-leg forms for all equations of the lists 
Q and H (the results for the list A follow from these ones). As it turns out, 
in almost all cases it is more convenient to write the three-leg equation (|35|) 
in the multiplicative form 

^(x,M;a)/^(a;, w;/3) = ^{x,y;a,P). (36) 

For the list Q the functions \E' and $ corresponding to the "short" and to 
the "long" legs, respectively, essentially coincide. One has in these cases: 

^{x,u;a) = F{X,U;A), ^{x,y; a, (3) = F{X,Y; A - B), (37) 

where some suitable point transformations of the field variables and of the 
parameters are introduced: x = f{X), u = f{U), y = f(Y), and a = p{A), 
/3 = piB). 

Theorem 7. The three-leg forms exist for all equations from Theorem 1. 
For the lists Q, H they are listed below. 

(Q1)5=q: An additive three-leg form with (j){x, y; a, [3) = ip{x, y;a — [3), 

ilj(x,u]a) = . (38) 

X — u 

{Ql)s=i: A multiplicative three- leg form with ^{x,y; a, P) = \E'(x, a— /?), 

^ix,u;a) = ^~^ + " . (39) 
X — u — a 

(Q2): A multiplicative three-leg form with $(x, y; a, (3) = \Ef(x, y^a — (3), 

T/ A {X + U + a){X ~U + a) 

^^^^^-^^^= \x + U-a)\x-U-ay ^''^ 
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where x = X"^, u = U'^ . 

{Q,3)s=o'- A multiplicative three- leg form with ^{x,y; a, (3) = '^{x,y\a/ (3), 

^{x,u]a) = — — -. (41) 
X — au 

Under the point transformations x = exp(2X), a = exp{2A), etc., there 
holds with 

(Q3)5=i: A multiplicative three-leg form with ^3l{), 

, sinhiX + U + A) sinUX - U + A) , , 
^[x,u;a) = F[X,\J;A)= . , ■ ^ — FT^-TT, 43 

^ ^ V ' J y smh(X - A) sinh(X -V - A) ^ ' 

where x = cosh(2X), a = exp{2A), etc. 

(Q4): A multiplicative three- leg form with j^Slj), 



^ r^r^ rr A^ (r{X + U + A)(t{X - U + A) 

-^(x,u:a) = F(XM:A) = —: ^-^ (44) 

^ ' ' ^ V , , ; a{X + U - A)a{X - U ~ A)' ^ ^ 

where x = p{X), a = p{A), etc. 
(HI): An additive three-leg form 

a- P 

^{x,u; a) = X -\- u, (f){x,y; a, P) = . (45) 

x-y 

(H2): A multiplicative three-leg form 

X — y-\-a — P 

'i/ix,u; a) = X -\- u -\- a, ^(x,y]a,pj = -. (46) 

X — y — a -\- p 

(H3): A multiplicative three-leg form 

^(x, u;a) = xu + 5a = exp{2X + 2U) + Sa, (47) 

^/ m px-ay sinh(X - f/ - A + 5) 

^{x,y-a,P) = ^ / = ■ , ^, 48 

ax — Py smn(A — U + A — B) 

where x = exp(2X), a = exp(2y4), etc. 

Proof. We start with the hst H, for which the situation is somewhat 
simpler. Finding the three-leg forms of Eqs. (HI) and (H3)^^q is almost 
immediate: these equations are equivalent to 

a-p u px-ay 

u — V = and — = --, (49) 

X — y V ax — Py 

respectively. For other equations of the list H one uses the following simple 
formula which will be also quoted on several occasions later on. 
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Lemma 8. The relation Q{x, u, v, y; a, /3) = yields 

h{x,u]a) 

h{x,v;p) Qu' 

Proof. 

h{x,u]a) _ g{x,u;a,P) _ QQy^ - QyQ^ _ 
h{x,v;P) g{x,v;P,a) QQyu - QyQu Qu 



(50) 



□ 



This lemma immediately yields the three-leg forms for the cases when 
Quv = 0: then the right-hand side of (pO|) does not depend on m, t> and can 
therefore be taken as y; a, while "^{x, u; a) := h{x, u; a). This covers 
the cases (H2) and (H3). Observe also that under the condition Quv = 
we have G{x,y; a, P) = QQuv — QuQv = —QuQv Thus, in these cases the 
function \1/ associated to the "short" legs is nothing but the polynomial /i, 
while the function $ associated to the "long" legs, is just a ratio of two linear 
factors of the quadratic polynomial G. 

The situation with the list Q is a bit more tricky. The inspection of for- 
mulas (ql), (q2), (q3), (q4) shows that in these cases h{x,u] a) is a quadratic 
polynomial in u, and, similarly, the polynomial G{x,y; a, (3) is a quadratic 
polynomial in y. We try a linear-fractional ansatz for $: 

^{x, u; a) = — -, $ X, y; a, (3) = — — , 51 

P-[x,u;a) s^[x,y;a,p) 

where the functions p±{x, u, a) are linear in u, and the functions s±{x, y; a, (3) 
are linear in and 

m; m; a) = h{x,u;a), 

s+{x,y;a,l3)s-{x,y;a,f3) = G{x,y; a, (3). 



According to (|2T1) , the coefficients of polynomials p±{-,u]-) and s±{-,y] , 
are rational functions of x and ^yr(x). This justifies "uniformizing" changes 
of variables, namely 

X = X"^ in the case (Q2), when r{x) = x, 

X = cosh(2X) in the case (Q3) with 6 = 1, when r{x) = — 1, 

X = p{X) in the case (Q4), when r{x) = 4x^ — g2X — g^. 

The ansatz ( [^ ) turns out to work. For the cases (Q1)-(Q3) identify the 
functions p± from (|5lD with the numerators and denominators of the expres- 
sions listed in (for the case (Q3)^^q take the fraction (0), i.e. the 
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expression through x, u). In the case (Q4) the functions p± hnear in u are ob- 
tained by dividing the numerators and denominators of the fractions in (|4^) 
by cr^(X)(T^([/). Make similar identifications for the functions s± from (^T]). 
So, under the point transformations used in the formulation of Theorem ^ 
we have: 



u; a) = P{X, U- ±A), s^{x, y- a, (3) = P{X, U; ±AtB), 



where 



ms=i 
m 
ms=o 
ms=i 

m 



P{x, u; a) = X — u + a, 
P(X, U;a) = {X + U + a){X - U + a), 
P(X, U; A) = sinh(X ~U + A), 
P{X, U; A) = sinh(X + U + A) sinh(X -U + A), 
a{X + U + A)a{X -U + A) 



PiX,U;A) 



a\X)a\U) 



A straightforward computation shows that 

p+{x, u; a)p-{x, v; y; a, (3) - u; a)p+{x, v; (3)s+{x, y; a, (3) 

= p{x;a,l3)Q{x,u,v,y;a,l3), (52) 

with some factor p depending only on x. (Obviously, the left-hand side 
of the latter equation is linear in u,v,y.) Concretely, we find: p = 2 in 
the case (Ql)5=i; p = iX in the case (Q2); p = x in the case (Q3)5=o; 
p = 2a/9sinh(2X) in the case (Q3)5=i; finally. 



p = <T\A)a\B) 



a{B - A) a{2X) 
g{B + A) ' a^{X) 



in the case (Q4). So, in all these cases Q = is equivalent to vanishing of 
the left-hand side of (|^), which, in turn, is equivalent to the multiplicative 
three-leg formula (^6]) with the ansatz (0). Finally, we comment on the 
origin of the additive three-leg form of (Ql)5=o- Rescaling in (|39| ) the pa- 
rameters as a I— 5a, we come to the three-leg form of the equation (Ql) 
with a general 5 7^ 0. Sending 5 ^ 0, we find: \Ef = 1 25iIj + 0(5^), 
$ = 1 + 250 + 0(5^), with the functions ip, cj) from (0), and thus we arrive 
at the additive formula for the case 5 = 0. □ 



Remark. The notion of the three-leg equation was formulated in [§. The 
three-leg forms for Eqs. (Ql), (Q3)5^q, (HI), (H3)^^q were also found there. 
The results for (Q2), (Q3)5^^, (Q4), (H2), (H3)^^j^ are given here for the first 
time. 
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6 Lagrangian structures 



Recall that the three-leg form of Eq. (|I]) imply that the discrete Toda 
type equation (0) holds on the black and on the white sublattices. (Of 



course, for multiplicative equations (|36D we set ijj{x,u;a) = log\l/(x,M;a 



(f){x,y; a, (3) = log a, /?).) We show now that all these Toda systems 

may be given a variational (Lagrangian) interpretation. Notice that only the 
functions 0(a;, a, /3) related to the "long" legs enter the Toda equations, 
and that the Toda systems for the cases (H1)-(H3) are the same as in the 
cases (Ql), (Q3)5=o. 

Consider the point changes of variables x = f{X) listed in Theorem 
in the cases (Ql), (HI), (H2), when no such substitution is listed, set just 
X = X. For the sake of notational simplicity, we write in the present section 
'?/'(x, m; a) for ip{f{X), f{U); a), etc. 

Lemma 9. For all equations from Theorem 1 there exist symmetric func- 
tions L{X, U; a) = L{U, X; a) and A{X, Y; a, (3) = A{Y, X; a, b) such that 

'ilj{x,u;a) = ij{f{X),f{U)-a) = ^L{X,U;a), (53) 
<P{x,y;a,P) = <P{fiX),fiYy,a,P) = -^A(X,Y; a, P). (54) 

Proof. It is sufficient to notice that the functions {d/dU)ip{x,u; a) are sym- 
metric with respect to the permutation X ^ U, and similarly for 0. □ 

This observation has the following immediate corollary. 

Proposition 10. Let be the "black" subgraph, and let E{^) be the set of 

its (non- oriented) edges. Let the pairs of labels {a, (3) be assigned to the edges 
from E{ir) according to Fig.^ so that, e.g., the pair (ai, 0:2) corresponds to 
the edge {x,Xi2). Then for all equations from Theorem 1 the discrete Toda 



type equations pl|j are the Euler-Lagrange equations for the action functional 



S= Yl A(X,F; «,/?). (55) 

{X,Y)eE{ir) 

This result could be anticipated, since it is quite natural to expect from 
a system consisting of equations on stars to have a variational origin. Our 
next result is, on the contrary, somewhat unexpected, since it gives a sort of a 
variational interpretation for the original system consisting of the equations 
on quadrilaterals (p. This is possible not on arbitrary quad-graphs but 
only on special regular lattices. We restrict ourselves here to the case of 
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the standard square lattice. For the cases (HI), (H3)^^q (the discrete KdV 
equation and the Hirota equation) our results coincide with those found in 
TU| , for all other equations from Theorem 1 they seem to be new. 

For the sake of convenience orient the square lattice as in Fig. ^ Denote 

Ml V2 

Pi. 




Figure 6: Two elementary quadrilaterals of the square lattice 



by El {E2) the subset of edges running from south-east to north-west (resp. 
from south-west to north-east). Add to the edges of the original quad-graph 
the set of the horizontal diagonals of all elementary quadrilaterals. Let 
the labels a be assigned to the edges from Ei, and the let the labels (3 be 
assigned to the edges from E2. It is natural to assign to the diagonals from 
E^ the pairs (a,/5) from the sides of the correspondent quadrilateral. 

Proposition 11. Solutions of all equations from Theorem 1 are critical for 
the following action functional: 

S= LiX,U;a)- L(X,V^;/?)- A(X,r;a,/3). (56) 

(X,U)€Ei {X,V)eE2 (X,y)G-B3 

Proof. For any vertex x, the correspondent Euler-Lagrange equation relates 
the vertices of two elementary quadrilaterals, as in Fig.|^. Due to Lemma |^, 
this equation reads: 

ijj{x, ui] ai) - 4j{x, vi] pi) - (f){x, yi; ai, + 
ip{x,U2;a2) -i'{x,V2;P2) - 4>{x,y2;a2,P2) = 0. 
This holds, since (^) is fulfilled on both elementary quadrilaterals. □ 

The variational interpretation allows one to find invariant symplectic struc- 
tures for reasonably posed Cauchy problems for Eqs. (|l]). As is well-known, 
one way to set the Cauchy problem is to prescribe the values x on the 
zigzag line like in Fig.|^ and to impose periodicity in the horizontal direc- 
tion (so that one is dealing with the square lattice on a cylinder) |T^, |TB| . 
Then Eq. (|lD defines the evolution in the vertical direction, i.e. the map 
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Figure 7: The Cauchy problem on a zigzag 



Proposition 12. Let the edges [xi, Xj+i) carry the labels ai, and let the edges 
(xjjJj+i) carry the labels cti, so thata2k = Ci2k-2, anda2k-i = C(2k+i- Denote 



d 92 
s{X, U;a) = Q^i^i^, u; a) = q^q^ L{X, U ; 



(57) 



so that s{X,U;a) = s{U, X; a) . Then the following relation holds for the 
map {xi} {xi}: 

s{Xi, Xi+i, ai)dXi A dXi+i = ^ s{Xi, X^+i; ai)dXi A rfXj+i . 



iG 



{5t 



Proof. We use the argument methodologically close to Consider the 

action functional § defined by the same formula as (|56D but with the sum- 
mations restricted to the edges depicted in Fig. and restricted to such fields 
X which satisfy the equation (|I]) on each elementary quadrilateral. Consider 
the differential 



dS 



all vertices x 



dX 



Due to Proposition |TT], the expression dS/dX vanishes for all vertices where 
six edges meet, i.e. for all vertices except X2k+i and X2k- So, we find: 



,53 \^ dL{X2k,X2k+i]a2k) dL{X2k+i,X2k+2]Oi2k+i, 
"S) — aA2fe+i 



k&LjNTL 



dX 



2k+l 



dX 



■dX. 



2k+l 



2k+l 



dL{X2k-l, X2k', 0!2k-l) ,:r^ dL{X2k, X2k+l', Ct2k) 

aA2fc H "A2fc. 



dX2k dX2k 
Differentiating this 1-form and taking into account (i^S = yields (|5B[). □ 
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Remarks. 1) Notice that only the functions related to the "short" legs 
are present in formula (|58|) . 

2) One is tempted to interpret formula ( |55D as symplecticity of the map 
{Xi] {Xj}. However, the 2-form ( ^81) is degenerate. One finds a genuine 
symplectic form, if one considers a quasi-periodic initial value problem in- 
stead of a periodic one, and extends the phase space by the correspondent 
monodromies (cf. |T^, for the case (H3)^^q (the Hirota equation)). Actu- 



ally, the formula (|58|) means the invariance of a degenerate 2-form only when 
Oii = Ui- This is easy to achieve by setting all aj = a, or, more generally, all 

tt2fc = Oi and all a2k+i = l3- 



We close this section with a local form of Proposition also based on 



the three-leg form of Eqs. (IT]). The particular case of the Hirota equation 



was given in [pT | 



Proposition 13. Associate the two-form uj{e) = s{X, U; a)dX/\dU to every 
(oriented) edge e = {x,u) of the quad-graph. Then Eq. ^ with the three-leg 
form ( Pi^ implies that the sum of uj 's along the boundary of any elementary 
quadrilateral, and therefore along any cycle homotopic to zero, vanishes. 



Proof. Differentiating (|35D and wedging the result with dX, we have: 

d d d 

-—i/j(x, u] a)dX AdU + -—i/j(x, v; (3)dV AdX = -—(pix, y; a)dX A dY. 
oU oV oY 

On the other hand, starting with the three-leg equation centered at y (or, in 
other words, flipping x ^ y, u v), we arrive at 

^^(y, v; a)dY A dV + ^V^(l/, u; (3)dU AdY = -^Hv, x; a)dY A dX. 

Adding these two equations, and taking into account Lemma ^ and the no- 
tation (^7D, we come to the statement of Proposition. □ 



The statement of Proposition is a particular case of Proposition 0, 
since the boundary of the domain in Fig.|^ (a cylindrical strip) is homotopic 
to zero. 



7 Relation to Backlund transformations 

In this section we interpret the integrable quad-graph equations of Theorem 
1 as nonlinear superposition principles (NSP) of Backlund transformations 
for the KdV-type equations. 
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First of all, we show how Backlund transformations themselves can be 
derived from our equations. Towards this aim, consider Eq. (|l]) on one 
vertical strip of the standard square lattice: 



S = {0,1} X Z. 

The fields on 5*0 = {0} x Z will be denoted by x^, while the fields on 5*1 = 
{1} X Z will be denoted by Uk- So, a single square of the strip S looks as 
on Fig.||. Suppose now that Xk = x{ke), /c e Z, where x{^) is a smooth 



Xk+l 



X Xf^ 



a 



a 



y = Uk+i 



U = Uk 



Figure 8: To the construction of Backlund transformations 



function, and similarly for Uk- In particular, this means that one has to set 
in Eq. (|l|) t> = x + ex^ + O(e^) and y = u + eu^ + O(e^). If, in addition, 
the parameter j3 is chosen properly, then Eq. approximates in the limit 
e ^ some differential equation which relates the functions x(^) and u{^). 
For the equations of the list Q the result is the most straightforward. 

Proposition 14. Set m Eq. (Ql) (3 = m Eq. (Q2) (3 = e'^/A, m Eq. 

(Q3) 13 = I- and in Eq. (Q4) (3 = p{e^). Then in the limit e ^ just 
described these equations tend to 

x^u^ = h{x,u; a), (59) 

with the correspondent polynomials h listed in the formulas (ql), (q2), (q3), 
(q4) of Proposition |^. 

Proof. The statement is almost obvious in the cases (Ql), (Q2). To see 
that it holds for Eq. (Q3), the latter can be first rewritten as 

/3(a^ — l)(x — v){u — y) = 

a{(3'^ - l){xv + uy) + (1 - (3){a'^ + (3){xy + uv) 
+{6^/Aa(3){a' - (3^){a^ - 1){(3^ - 1). 



30 



If /? = 1 — then the above equation approximates 

(a^ - l)x!:U!: = -a{x'^ + u^) + (a^ + l)xu - {6'^/Aa){a'^ - 1)^. 
Finally, in the most intricate case (Q4) one starts by rewriting the equation 



as 



— — — {x — v){u — y) = xuvy H — ^ {xuv + uvy + vyx + yxu) 

IQjQ CLq 

' ~'~ '^^ (x + f ) (m + y) + — {xv + uy) + — {x + u + V + y) + 



2,CLq CLq CLq CLq 

From expressions for the coefficients Oj given in Theorem it is easy to see 
that if /3 = p(e^) ~ e~^, so that b = p'{e^) ~ — 2e~^, then the left-hand side 
of the above equation tends to ax^u^, while the right-hand side tends to 

x'^u^ — 2axu{x + u) — ^2a^ — ^^xu + + v^) 

+ (g?, + y a) + ^) + (H + fl'3") = <^K^^ a). 

This proves the proposition. □ 
Eq. (0), read as a Riccati equation for u with the coefficients depen- 
dent on X, describes a transformation x ^ u, which turns out to be a 
Backlund transformation for the Krichever-Novikov equation ^0 : 



3 

xt = x^^^ - - (60) 

with the polynomial r{x) being the discriminant of h{x, u; a). In other words, 
if X is a solution of (^), and u is related to x by (^9|), then u is also a 
solution of (0). It should be noticed that, in turn, the partial differential 
equation (|60D may be derived from (^), either through a sort of continuous 
limit, or as a higher symmetry. In any way, it would be fair to say that the 
whole theory of the equation ( |60D and its Backlund transformations (|59D is 
contained in the correspondent quad-graph equation (|Tp from the list Q. To 
complete the picture, we demonstrate that all equations (|l]) listed in Theorem 
1, in turn, can be interpreted as nonlinear superposition principles for the 
Backlund transformations. To this end consider the system of four differential 
equations of the type (|59D corresponding to four sides of the quadrilateral on 

x^u^ = h{x,u; a), u^y^ = h{u,y; P), (61) 
x^v^ = h{x, v; /?), v^y^ = h{v, y; a). (62) 
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We will consider it for functions h corresponding to all cases listed in Theorem 
1, except for those two when h actually does not depend on parameters, 
namely (HI) with h{x,u) = 1, and (H3)^^q with h{x,u) = xu. In all other 
cases the system (|6lD, (|6^) makes perfect sense and its consistency is quite 
nontrivial. 

Proposition 15. The equation Q{x,u,v,y; a, f3) = is a sufficient condi- 
tion for the consistency of the system of differential equations ^\), ^6^). 
Moreover, it is compatible with this system, i.e. 

QdQ=o = 0. (63) 

Proof. The consistency condition of the differential equations (|6ll ) , (|6^) reads: 

Q := h{x, u; a)h{v, y; a) — h{x, v; P)h{u, y; (3) = 0. (64) 

The left-hand side Q of this equation is a polynomial of degree 1 in each 
variable for the equations of the list H, and of degree 2 in each variable 
for the equations of the list Q. In the former case it is directly seen that 



Eq. (]6^) exactly coincides with Q{x,u,v,y; a, /3) = 0. In the latter case 
it can be shown that the polynomial Q divides the left-hand side of (p^. 
More precisely, it is verified that, up to a constant factor, Q = QP, where 
in the cases (Ql), (Q2) one has P = (5|/3^-/3, in the case (Q3) one has 
P = Q\f3^i/f3, and in the case (Q4) one has P = Q|(/g^b)^(^^_b). This proves 
the first statement of proposition. The second one reads as 

{QxX^ + QuU( + QvV^ + Qyy^)\Q=o = 0. 

We will prove that actually both terms QxX^ + Qyy^^ and QuU^ + Q^f^ vanish 
separately. For example, 

{QuU^ + QvV^)\q=o = ^ {Quh{x,u;a) + Qyh{x,v; = 0. 

The last step follows from Lemma ^ □ 



In the two cases when the system (|6lD, ( |62D becomes trivial, there still 
exist Backlund transformations of a different kind, such that the equation 
Q = serves as their NSP. Namely, in the cases (HI), (H3)^^q the following 
equations come to replace (^): 

x^ + = (x — uY + a, (65) 

xc U( a /X u\ , , 

- + - = (- + -), 66 

X U 1 \U X/ 
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respectively. The second of these equations is probably better known in the 
coordinates x = exp(X), u = exp(f/): 

X^ + U^ = a cosh(X - U) . (67) 

Eq. ( |B3D defines a Backlund transformation for the potential KdV equation 

a:^ = a;^^^ _ ex^. (68) 

Similarly, Eq. ( |67|) defines a Backlund transformation for the potential 
MKdV equation 

Xt = X^^^ - 2X|, (69) 

or, alternatively, for the sinh-Gordon equation which belongs to the same 
hierarchy. 

8 Conclusions and perspectives. 

Three— dimensional consistency as integrability criterium; Yang- 
Baxter maps. One of the traditional but rather ad hoc definitions of the 
integrability of two-dimensional systems is based on the notion of the zero- 
curvature representation. For a system on a quad-graph consisting of equa- 
tions (|ip with fields associated to the vertices of the elemntary quadrilateral 
on Fig. |l|, the zero curvature representation is usually encoded in the formula 
like 

L{y, u, P; X)L{u, x, a; A) = L{y, v, a; A)L(t>, x, P; A), 

where A is the spectral parameter, so that the matrices L take values in 
some loop group. In and independently in it was demonstrated how 
to derive the zero curvature representation from the three-dimensional con- 
sistency. It should be mentioned, however, that to assign fields to the vertices 
is not the only possibility. Another large class of two-dimensional systems 
on quad-graphs build those with the fields assigned to the edges, see Fig. 
^. In this situation it is natural to assume that each elementary quadrilat- 
eral carries a map R : X'^ X'^, where X is the space where the fields 
a,b take values, so that (a2,&i) = R{a,b; a, P). The question on the three- 
dimensional consistency of such maps is also legitimate and, moreover, began 
to be studied recently. The corresponding property can be encoded in the 
formula 

^23 O ^13 O Rl2 = Rl2 O ^13 O ^23, (70) 
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Figure 9: An elementary Figure 10: Three-dimensional consis- 

quadrilateral; both fields and tency; fields assigned to edges 

labels are assigned to edges 



where each Rij : t— > acts as the map R on the factors z,j of the 
cartesian product X^ and acts identically on the third factor. This equation 
should be understood as follows. The fields a, h are supposed to be attached 
to the edges parallel to the 1st and the 2nd coordinate axes, respectively. 
Additionally, consider the fields c attached to the edges parallel to the 3rd 
coordinate axis. Then the left-hand side of ( [70| ) corresponds to the chain of 
maps along the three hind faces of the cube on Fig. |l^: 

(a, h) ^^ (a2, 61), (a2, c) (023, ci), (61, ci) ^^ (613, C12), 

while its right-hand side corresponds to the chain of the maps along the three 
front faces of the cube: 

(6, c) ^^ (63, C2), (a, C2) ^ (03, C12), (03, ^3) ^ (023, &13) 

So, Eq. (|70D assures that two ways of obtaining (023, &13, C12) from the initial 
data (a, 6, c) lead to the same results. The maps with this property were 



introduced by Drinfeld under the name of "set-theoretical solutions of 
the Yang-Baxter equation", an alternative name is "Yang-Baxter maps" 
used by Veselov in the recent study |3^, see also references therein. Under 



some circumstances, systems with the fields on vertices can be regarded as 
systems with the fields on edges or vice versa (this is the case, e.g., for the 
systems (Ql), (Q3)5=o, (HI), (H3)5=o of our list, for which the variables X 
enter only in combinations like X — U for edges (x, u)), but in general the two 
classes of systems should be considered as different. The notion of the zero 
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curvature representation makes perfect sense for Yang-Baxter maps: such a 
map can be called integrable, if it is equivalent to 



L{b, (3\ X)L{a, a; A) = L{c, a; X)L{d, (3; A). 
The problem of integrability of Yang-Baxter maps in the sense of existence 



of a zero-curvature representation is under current investigation [^. Also 
the problem of classification of Yang-Baxter maps, like the one achieved in 
the present paper, is of a great importance and interest. 

A different direction for the development of the ideas of the present 
paper constitute quantum systems, or, more generally, systems with non- 
commutative variables. To remain in the frame of the present paper, these 
are systems (0), where the fields {x,u,v,y) take values in an arbitrary asso- 
ciative (not necessary commutative) algebra with a unit. It turns out that 
the notion of the three-dimensional consistency can be formulated also for 
such non-commutative systems. Also the derivation of the zero curvature 
representation can be extended to the non-commutative framework 

It should be mentioned that in the area of the three-dimensional consis- 
tency of classical systems there also remains a number of interesting open 
problems. For instance, one of the assumptions under which the classification 
was carried out in the present paper, was less natural, namely the tetrahe- 
dron condition. As we pointed out in Sect. |^, there exist three-dimensionally 
consistent equations without the tetrahedron property, however all examples 
we are aware of are trivial (linear or linearizable) : 



Q{x, u,v,y) = x + y — u — v = or Q{x, u, v,y) = xy — uv = 



or those obtained from these two by the action of a Mobius transformation 
on all variables. These examples do not contain parameters, and thus the 
three-dimensional consistency does not give a zero curvature representation 
with a spectral parameter for them (their integrability is anyway obvious). 
It would be interesting to find out whether there exist nontrivial examples 
violating the tetrahedron property. 

There is also a vast field of multi-field integrable equations on quad- 
graphs. Existing examples indicate that their study is very promising. 

Four— dimensional consistency of three— dimensional systems. Ve- 
ry promising is also the application of the consistency approach to the three- 
dimensional integrability. The role of an ad hoc definition of integrability, 
played in two dimensions by the zero curvature representation, now goes to 
the so called local Yang-Baxter equation, introduced by Maillet and Nijhoff 
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22| . The role of the transition matrices from the zero curvature representa- 
tion is played in this novel structure by certain tensors attached to the ele- 
mentary two-dimensional plaquettes of the three-dimensional lattice. There 
exist a number of results on finding this sort of structure for some three- 



dimensional integrable systems [ll9| , {Wj , [17| . It would be desirable to relate 



this ad hoc notion of integrability to some constructive one. In the spirit of 
the present paper, this constructive notion should be the four-dimensional 
consistency. 

In the three-dimensional context there are a priori many kinds of systems, 
according to where the fields are defined: on the vertices, on the edges, or 
on the elementary squares of the cubic lattice. Consider first the situation 
when the fields are sitting on the the elementary squares. Attach the fields 
a, b, c to the two-dimensional faces parallel to the coordinate planes 12, 23, 
13, respectively, so that a, b, c are sitting on the bottom, the left and the 
front faces of a cube Fig. |^, and 03, bi, 02 on the top, the right and the back 
faces. The system under consideration is a map S : ^ attached to the 
cube, so that S{a,b,c) = (03,61,02). The condition of the four-dimensional 
consistency of such a map can be encoded in the formula 

5*134 O 5*234 O 5*124 ° 5*123 = '5*123 ° '5*124 ° '5*234 ° 5*i34. (71) 

This equation should be understood as follows. Additionally to the fields a, b, 
c, consider the fields d, e, /, attached to the two-dimensional faces parallel to 
the coordinate planes 24, 14, 34 of the four-dimensional hypercubic lattice. 
Each map of the type Sijk in ( |7TD is a map on X^{a,b,c,d,e, f ) acting as 
5 on the factors of the cartesian product corresponding to the variables 
sitting on the faces parallel to the planes ij, jk, ik, and acting trivially on 
the other three factors. Thus the left-hand side of ( [71]) corresponds to the 
chain of maps 

(a, b, c) ^^ {as, bi, C2), (03, d, e) ^ (034, di, 62), 

(61, rfi, /) (614, rfi3, /2), (C2, 62, /2) ^ (C24, 623, /21), 

while the right-hand side of (0) corresponds to the chain of maps 

(c, e, /) ^-> (c4, 63, /i), (6, d, fi) ^ (64, ds, /12), 
(a, 4, 63) ^^ (04, di3, 623), (04, 64, C4) ^ (034, bu, C24). 

Eq. (^1]) expresses then the fact that two different ways of obtaining the data 
(034, 614, C24, (ii3, 623, /12) from the initial data (a, 6, e, d, 6, /) lead to identical 
results. This equation is known in the literature as the functional tetrahedron 
equation |T^, . (Note that the standard notation used in the literature on 
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the tetrahedron equation is different: the indices 1 < 7 < 6 of Sap-^ 



numerate the two-dimensional coordinate planes.) The paper contains 
also a list of solutions of this equation with A" = C, possessing local Yang- 
Baxter representations with a certain ansatz for the participating tensors. 
One of the most remarkable examples is the star-triangle map: 

a , b c 

"3 = 1 1 ' ^1 = , C2 = — - — . 72 

ab — be — ca ab — be — ea ab — be — ea 

(Usually this equation is written in a more symmetric form which is obtained 
by changing c t— > — c, with all plus signs in all denominators; however in that 



form the map does not satisfy Eq. ([711).) See also |T^ for an alternative 
local Yang-Baxter representation for this system. In |3l|, ^ an important 
solution of the functional tetrahedron equation with A* = was introduced 
and studied in detail. There seems to be no method available for deriving the 
local Yang-Baxter representation for a given map satisfying the functional 
tetrahedron equation. 

Further, consider three-dimensional systems with the fields sitting on the 
vertices. In this case each elementary cube carries just one equation 

Q{x, xi, X2, X3, X12, X23, X123) = 0, (73) 

relating the fields in all its vertices. Such an equation should be solvable for 
any of its arguments in terms of the other seven ones. The four-dimensional 
consistency of such equations is defined as follows. 

• Starting with the initial data x, Xi (1 < i < 4), (1 < i < j < 
4), the equation ( [73| ) allows us to uniquely determine all fields Xijk 
{l<i<j<k<A). Then we have four different possibilities to find 
2^1234, corresponding to four three-dimensional cubic faces adjacent to 
the vertex 0:1234 of the four-dimensional hypercube, see Fig. p. All four 
values actually coincide. 

So, one can consistently impose equations (|7^) on all three-dimensional cubes 
of the lattice Z^. It is tempting to accept the four-dimensional consistency of 
equations of the type ( [73| ) as the constructive definition of their integrability. 
It will be very important to solve the correspondent classification problem. 
We expect that this definition will allow one to derive the local Yang-Baxter 
representation, as a replacement for the zero curvature representation char- 
acteristic for the two-dimensional integrability. 

We give here some examples. Consider the equation 

(Xi - X3){X2 - X123) _ (X - Xi3)(xi2 - X23) 
(X3 - X2)(xi23 - Xi) (Xi3 - Xi2)(x23 - x) ' 
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This equation appeared for the first time in |29|, |T8|, along with a geometric 



interpretation. It is not difficuh to see that Eq. (jTj) admits a symmetry 
group of the cube. This equation can be uniquely solved for a field at 
an arbitrary vertex of a three-dimensional cube, provided the fields at other 
seven vertices are known. The fundamental fact not mentioned in |jl8| is: 

• Eq. (|7^) is four-dimensionally consistent in the above sense. 

This is closely related (in fact, almost synonymous) to the functional tetra- 
hedron equation for the star-triangle map ([7^), since the plaquette variables 
(a, 6, c) of the latter can be "factorized" into combinations of vertex variables 
X of Eq. (|73p. More precisely, given a solution of (^) and setting 

X12-X X23 -X Xi'i- X 

a = , = , c = 

Xi - X2 X2- X3 Xi - Xs 

we arrive at a solution of (^). 

A different "factorization" of the plaquette variables into the vertex ones 
leads to another remarkable three-dimensional system known as the discrete 
BKP equation [^, |jl8|. For any solution a; : 1— >■ C of ([74|), define a 



function r : Z'^ t— > C by the equations 



TjTj Xij X . /^7^^ 

— 1 < 3- (75) 



TT^j Xi X') 



1 



Eq. (0) assures that this can be done in an essentially unique way (up 
to initial data on coordinate axes whose influence is a trivial scaling of the 
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solution). The function r satisfies on any three-dimensional cube the discrete 
BKP equation: 

TTijk - TiTjk + TjTik - TkTij = 0, i < j < k. (76) 

The following holds: 

• Eq. (|76D is four-dimensionally consistent. Moreover, for the value T1234 
one finds a remarkable equation: 

Tri234 - 7-12^-34 + T13T24 - r23r34 = 0, (77) 

which essentially reproduces the discrete BKP equation. So, ri234 does 
not actually depend on the values Tj, 1 < i < 4. This can be considered 
as an analog of the tetrahedron property of Sect. 0. 

Notice that usually the discrete BKP equation (^) is written in a slightly 
different and more symmetric form, with all plus signs on the left-hand side. 
On every three-dimensional subspace these two forms are easily transformed 
into one another. However, this cannot be done on the whole of Z^. Eq. ([76| ) 
with all plus signs on the left-hand side does not possess the property of the 
four-dimensional consistency. 

Further, we mention systems of the geometrical origin (discrete analogs of 
conjugate and orthogonal coordinate systems) [|n|], |^, which also have the 
property of the four-dimensional consistency. We plan to address various 
aspects of the four-dimensional consistency in our future publications. 
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